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Preceding 


INTRODUCTION 


This  study  was  undertaken  at  the  request  of  the  Selected  Systems  Effec¬ 
tiveness  Program  (SSEP)  of  the  Joinl  Technical  Coordinating  Group  for 
Munitions  Effectiveness  (JTCG/ME).  The  study  is  an  effort  to  determine 
the  effectiveness  of  present  Department  of  Defense  methodology  in  de¬ 
termining  or  predicting  the  diffusion  of  chemicals  and/or  radiological 
matter  in  the  atmosphere.  In  par-icula'*,  the  study  shows  the  extent 
to  which  a  difference  in  temperaiure  between  two  levels  i car  the  ground 
(static  stability)  or  dynamic  stability  (ratio  z/L)  is  r-elated  to  the 
govern i..g  parameters  of  turbulence  in  the  diffusion  models  considered. 


DISCUSSION 


Haugen  and  F.jquay  C  j  pointed  out  that  one  should  consider  several 
meteorological  parameters  in  characterizing  the  diffusi/e  power  of  the 
atmosphere.  Among  these  parameters  are  the  following: 

I.  The  mean  wind  velocity,  determining  trie  path  of  the  cloud  jown- 
wind  from  the  source  and  the  distribution  of  particles  along  the  mean 
wind  path.  In  the  case  of  a  continuous  source  this  downwind  transport 
results  in  a  lower  concentration  of  particles  as  the  mean  wind  increases. 


2.  The  RMS  (standard  deviation)  of  The  wind  direction  determin¬ 

ing  the  shape  and  distortion  of  the  cloud,  and  in  a  sense  serving  as  an 
indication  of  horizontal  mixing. 


t  •  r»»  ii'  <  .  i  ,  .1 «  .  j  j  _  •  \  x_l  —  .  .  x  t  | 
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( Ojr ) ,  determining  the  dispersion  of  the  cloud  in  the  vertical. 


4.  The  vertical  temperature  gradient  (AT)  as  an  indicator  of  verti¬ 
cal  mixing.  If  the  temperature  increases  with  height,  the  vertical  mixing 
is  limited  to  mechur ical  turbulence  transport,  whereas  both  mechanical 
and  convective  mixing  generally  occur  during  nen-inversion  conditions. 


Relationships  between  the  four  parameters  discussed  above  have  been 
examined  by  Haugen  and  Fuquay  Clj  and  Record  et  al.  [2j.  Models  which 
quanti fatively  predict  Turbulent  diffusion  often  assume  that  these 
relationships  exist.  In  particular,  the  United  States  Army  Combat  Devel¬ 
opment  Command's  Field  Manual  3-10,  "Employment  of  Chemical  and  Biological 
Agents,"  utilizes  a  AT  between  heights  of  0.5  and  4.0  meters  to  determine 
the  horizontal  and  vertical  rates  of  iurbulent  mixing.  Various  levels 
have  been  used  for  the  measurement  of  the  temperature  gradient;  e.g., 
project  Ocean  Breeze  and  Dry  Gulch  [JD  used  temperatures  at  heights  of 
2  and  17  meters.  More  recently.  Record  et  al.  [2]  used  heights  of 
approximately  3  and  20  meters.  The  actual  levels  at  which  these  tempera¬ 
ture  measurements  were  made  were  based  primarily  on  existing  instrumentation 
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at  the  site.  Because  of  the  highly  varying  conditions  observed  near 
the  surface,  there  is  a  need  to  investigate  the  effects  of  choice  of 
heights  on  turbulent  diffusion  parameters. 

SOME  FUNDAMENTAL  CONCEPTS 

A.  Wind  and  Its  Variability 

Wind  is  i  "Hated  by  pressure  forces  which  are  created  by  density  gradi¬ 
ents  resulting  from  differential  heating.  I:  \z  "T  d  i  H  ed  by  tne  rotation 
of  the  earth,  frictional  effects,  and  bv  centrifugal  forces  due  tc  +he 
curvature  of  the  air  trajectory.  In  the  surface  boundary  layer  curvature 
and  coriolis  effects  often  may  be  neglected.  The  speed  and  direction 
of  the  wind  then  depends  primarily  on  the  magnitude  and  direction  of  the 
pressure  gradient  force,  and  on  the  properties  of  the  surface  which 
determine  the  magnitude  of  the  frictional  force. 

The  influence  of  friction  decreases  with  altitude  and  may  vary  with 
horizontal  distance.  When  the  surface  is  rough,  the  wind  speed  increases 
rapidly  at  heights  just  a  few  meters  above  the  ground  and  less  rapialy 
at  greater  heights.  When  the  ground  is  smooth,  the  increase  of  wind 
speed  with  height  is  less  pronounced  but  is  still  greater  in  the  lowest 
few  meters.  Mechanical  turbulence  created  by  roughness  elements  decreases 
with  height  more  rapidly  over  a  rough  surface  than  a  smoo+h  surface. 

When  the  degree  of  surface  roughness  changes  horizontally,  changes  in 
wind  and  the  intensify  of  turpulence  at  a  given  height  result.  Surface 
features  such  as  hillocks,  vegetation,  and  changes  in  other  lerrain  fea¬ 
tures  such  as  soil  composition  (which  may  lead  to  differential  heatinq 
and  variations  in  the  pressure  gradient  force)  Induce  spatial  variations 
in  wind  speed  and  turbulence. 

At  present,  it  is  not  possible  to  determine  analytically  the  micro-  and 
mesoscale  forces  in  space  and  time  which  affect  the  wind.  It  is  possible, 
however,  to  determine  the  average  magnitude  of  the  forces  in  space,  and 
hence  the  average  wind  may  be  determined  with  reasonable  accuracy.  Varia¬ 
tions  in  the  wind  in  both  space  and  time  within  local  areas  are  usually 
determined  statistically,  leaving  much  to  be  desired  in  terms  of  accuracy. 
It  is  this  variability  in  wind  which  leads  to  extreme  complexities  in 
problems  of  diffusion  and  transport  of  airborne  substances,  impact  pre¬ 
dictions  of  unguided  rockets,  and  numerous  and  varied  meteorological 
problems,  sucn  as  heat  and  momentum  transfer  and  the  coupling  actions 
between  layers  of  the  atmosphere  with  differing  wind  speeds  and  stabiliiy. 

B.  Definition  of  Turbulence 

The  vector  wind  can  be  written  as 

V  =  V  +  v'  (I) 
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where  the  bar  denotes  an  average  value,  an  arrow  3  vector  quantity  (the 
absence  of  an  arrow  denotes  a  scalar  quantity),  and  the  prime  a  devia¬ 
tion  from  the  average.  Equation  (I)  may  be  written  in  component  form  as 


u  =  u  +  u' 


v  =  v  +  v’ 


(2) 


w  =  w  +  w’ 


where  u,  v,  and  w  refer  to  the  components  of  the  wind  along  the  ortho¬ 
gonal  axes  x,  y,  and  z,  respectively.  Deviations  from  the  average  are 
referred  to  is  turbulence  and  represent  that  portion  of  the  wind  which 
must  be  treared  staristical ly.  The  average  may  or  may  not  be  a  func¬ 
tion  of  time.  A  hypothetical  time  trace  of  wind  speed  showing  u  and  u* 
at  a  given  location  near  the  ground  is  illustrated  in  Figure  i.  The 
variations  in  spted  are  caused  by  variations  in  the  forces  discussed 
above. 

C.  Factors  ResDonsible  for  Turbulence 


A  convenient  way  of  looking  at  the  rate  of  growth  and  decay  of  turbulent 
energy  is  tc  examine  the  terms  in  the  eddy  kinetic  energy  equation. 
Lumley  and  Panofsky  C3H  wrote  this  equation  in  the  simplified  form 


(A)  (B) 


D(KE')  =  pK 
Dt 


(C) 

Kh  pg  36 
0*  3z 


(D) 

-  e 


(3) 


Here  (KiT' )  is  the  average  kinetic  energy  of  the  three  components  of 
turbulence,  Kj ,  and  ana  the  eddy  exchange  coefficients  for  heaT  and 
momentum,  respectively,  6  is  an  average  potential  temperature  in  degrees 
Kelvin,  p"  is  an  average  density,  g  is  gravity,  ^  is  the  mean  wind  velo¬ 
city,  z  is  height,  and e  represents  the  dissipation  cf  mechanical  energy 
into  heat.  The  terms  comprising  this  equation  are  interpreted  as  f o' lows 

I.  Term  (A)  represents  the  time-rate-of -change  of  the  average 
turbulent  energy  following  the  mean  motion.  Under  equilibrium  conditions 
and  neglecting  any  advective  processes,  this  term  will  be  zero.  A  posi¬ 
tive  value  of  term  (A)  represents  the  growth  and  a  negative  value  the 
decay  of  turbulent  energy. 
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2.  Term  (B)  represents  the  work  done  by  the  turoulent  stresses 
against  the  rates  of  mean  strain.  This  term  is  usually  positive  and 
thus  contributes  to  the  growth  of  turbulent  energy. 

3.  Term  (C)  represents  the  work  done  against  buoyancy.  This  term 
can  be  an  energy  source  or  sink  depending  on  the  gradient  of  potential 
temperature,  and  it  may  be  negligible  when  the  gradient  vanishes. 

4.  Term  (D)  represents  the  rate  of  dissipation  of  turbulent  kinetic 
erergy.  Since  dissipation  is  always  a  positive  quantity,  the  negative 
sign  makes  the  onl i re  term  an  energy  sink. 

Richardson  (see  Sutton  E4])  postulated  that  fhe  state  of  turbulence  can 
be  determined  from  the  ratio  of  term  C  to  term  B  in  Eq.  (3),  viz.. 


Rf  =  9  39/9z _  (4) 

^  9  ( | 3V/3z | )2 

Rf  is  referred  to  as  the  flux  Richardson  number.  If  it  is  assumed  that 
K^/l^  =  I,  then  Eq.  (4)  becomes  the  gradient  Richardson  number,  Ri.  If 
Rf  =  I ,  then  the  buoyancy  forces  remove  turbulent  energy  at  the  same 
fqte  that  it  is  being  produced  by  the  shearing  stresses.  This  does  not 
mean  i.>ai  Rf  =  I  is  a  critical  value  for  predicting  the  onset  of  turbulent 
moticn  but  rather  the  condition  for  equilibrium.  It  is  clear,  however, 
that  Rf  <  I  indicates  turbulent  energy  growth  and  Rf  >  I  indicates  decay. 

If  is  also  possible  to  express  Rf  as  a  function  of  height  fn  non- 
dimensional  form  using  similarity  theory  C3H.  A  non-dimensional  wind 
shear,  S,  can  be  defined  as 


c  kz  3V  (5) 

s  =  W  dl 

where  k  is  von  Karman's  constant,  u*  is  the  fricT’on  velocity,  and  the 
other  variables  are  as  previously  defined.  A  scaiing  length,  L,  defined 
by  Mon  in  and  Obukhov  (see  Lumlcy  and  Panofsky  D]),  is  given  by 


CDp0u*3 

L  5  -  _  (6) 

kgH 

where  Cp  is  the  specific  heat  of  air  at  constant  pressure,  and  H  is  the 
vertical  flux  of  heat.  Using  Eqs.  (4),  (5),  and  (6)  and  tne  following 
expressions  for  H  and 
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hh-pcpkh  || 


(7) 


Xm 


u1 


*2/3v_ 
3z 


V- 


we  arrive  at 


SRf  =  _ =  z/L  (8) 

CppSu*3 

Observation  has  shown  C3H  that  in  near-neutrai  conditions  S  =  I  so  that 
Rf  =  z/L.  Hence,  z/L  is  a  measure  of  dynamic  stability  and  thereby  is 
a  means  of  ascertaining  whether  or  not  there  is  a  growth  or  decay  of 
turbulent  kinetic  energy. 

D.  Diffusion  Theory 

The  intensity  of  turbulence,  I,  is  given  by  Slade  [5D  as 


where  the  a's  represent  the  standard  deviation  of  the  longitudinal  (u), 
lateral  (v),  and  vertical  (w)  components  of  turbulence,  V  is  the  average 
wind  speed,  and  x,  y,  and  z  refer  to  coordinate  directions.  Longitudinal 
refers  to  the  direction  along  the  mean  wind,  and  lateral  to  the  per¬ 
pendicular  to  the  mean  wind  direction.  The  diffusing  power  of  the  atmo¬ 
sphere  is  directly  related  to  ttie  I ' s.  When  the  average  wnd  direction 
is  along  the  x-axis  such  that  v  =  0,  then  there  is  no  y  component  and 
we  can  write 


°A  = 


and 


(10) 


where  A  is  wind  direction  (azimuth)  in  radians,  and  a<5>  is  the  standard 
deviation  of  the  inclination  of  the  vector  wind  in  radians  (when  0=0, 
the  wind  is  horizontal).  In  practical  diffusion  work  the  basic  problem 
is  to  relate  the  measurable  quantities  given  in  Eqs.  (9)  and  (10)  to 
the  dispersion  of  airborne  effluents  under  various  stability  conditions. 
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Statistical  +heories  of  turouient  diffusion  show  that  the  variance  of 
particle  diffusion  in  the  y-direction  under  different  condito.is  is 
given  by  f63 


y2(-f;  =  2Kt  (t  large) 


(i  I) 


2 

y 


(t) 


(t  sma  1  i ) 


(12) 


and 


2(+)  =  1/2  Cy  (ut)2_n  (t  intermediate)  (13) 

where  Gaussian  distribution  of  the  particles  has  been  assumed,  and  where 
y(t)  =  cross  wind  distance  that  a  particle  moves 
from  some  origin 
K  -  eddy  diffusion  coefficient 
t  =  time 

Cy  =  eddy  diffusion  parameter  in  the  y-direction 
n  =  stability  parameter 
v'  -  crosswind  turbulent  fluctuation 
u  =  mean  wind  speed 

The  constant  K  in  Eq.  (II)  and  Cy  in  Eq.  (>3)  depend  upon  the  intensity 
of  turbulence  as  given  in  Eqs.  (S)  and  (iO).  While  it  will  not  be  shown 


the  variance  of  the  lateral  velocity  fluctuations,  and  stability.  The 
influence  of  stability  is  contained  in  the  parameter  n  defined  from  the 
wind  profile  by  the  relation 
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_  =  I  I 


(2-n) 


v2 


r^l 


(14) 


where  Vj  is  measured  at  z.  and  V2  at  22*  Different  vajues  of  n  represent 
different  stability  conditions.  Thus,  the  value  of  y2  in  Eqs.  (II)  -  (13), 
which  determines  the  distribution  of  particles,  in  the  cross-wind  direction 
(along  the  y-axis),  is  related  to  .the  intensity  of  turbulence  given  by  Eq. 
(9),  the  total  energy  spectrum  of  turbulence,  variability  of  wind  direc-  • 
tion  (Eq.  (10)),  stability,  and  average  wind  speed.  As  discussed  pre¬ 
viously,  these  variables  are  interrelated  with  each  one  usually  expressed 
as  a  function  of  stability.  Thus,  the  rate  of  diffusion  is  related  to 
the  degree  of  stability. 


E.  Meteorological  Variables  Related  to  Diffusion  Parameters 

The  meteorological  variables  most  commonly  related  to  atmospheric  diffu-  1 
sion  near  the  ground  include:  i)  tne  intensity  of  turbulence,  2)  the 
spectrum  of  turbulence,'  3)  variability  of  wind  direction  (vertical  and 
horizontal),  4)  stability,  and  5)  wind  speed.  These  parameters  are  re¬ 
lated  in  some  way  to  K  and  Cy  in  Eqs.  (IP  and  (13),  or  to  the  produc¬ 
tion  terms  in  the  kinetic  energy  equation  for  turbulence  (.Eq.  (3)). 

Static  stability  and  wind  speed  near  the, ground  are  related  through  Eq* 
(14).  To  illustrate  the  relatjonship  between  stability  and  the  diffusion 
coefficients,  the  following  data  were  taken  from  page  243  of  Atmospheric 
Diffusion  by  Pasqui I  I  [73. 


TABLE 


DIFFUSION  PARAMETERS  VS.  ATMOSPHERIC  STABILITY 


Stabi 1 i tv 

t  -  / °c \ 

1 18?  ft  *5  ft  '  ' 

Pi 

C  (mn/2) 

y 

Neutra 1 

‘  -1  and  0 

0.25 

0.095-0. 14 

Moderate  inversion 

+1  to  +3 

0.55 

0.052-0.077  , 

Strong  1 nvers i on 

✓*■ 

O 

0.50 

0.029-0.074 
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lhe  diffusron  coefficient  in  the  lateral  direction,  Cw,  is  assumed 
to  equal  that  in  the  vertical  direction.  Static  stability  is  a  function 
of  wind  spped 9  tsr«r2 i n  cond i o,n,s ?  "fho  *fypo  of  sir  moss  prosorsf’^  cloudi** 
ness  (or  radiation),  and  other  parameters.  It  is  observed  that  the  con¬ 
centration  of  pollutants  increases  over  cities  during  stable  conditions 
v:h?le  rapid  dispersion  takes  place  during  unstable  conditions.  While  the 
magnitude  of  the  diffusion  coefficient  is  proportional  to  instability, 
the  relationship  between  the  concentration  of  pollutants  and  stability 
conditions  does  not  always  exist.  For  example,  mechanical  turbu lance  may 
exist  even  during  inversion  conditions,  which  leads  to  an  increase  in  Cy. 

A  wide  range  of  values  for  the  diffusion  coefficients  under  different 
atmospheric  conditions  are  summarized  by  Pasqnill  Tl~\  and  Gifford  T6T. 

A  rather  wide  range  of  values  has  been  obtained  from  experimental  results. 


DA  l  A 

The  data  utilized  in  this  study  were  obtained  from  instruments  set  up  in 
an  array  forming  a  "T"  over  distances  of  300  meters.  The  instruments 
ut i  I  i zed  were  R.  m.  iouny  company^  UVW  (propel  ler  anemometers)  along 
with  unaspirated  thin  wire  thermocouples.  The  accuracies  of  these  instru¬ 
ments  and  data,  spacing  of  poles  on  which  instruments  were  mounted,  arid 
terrain  characteristics  have  been  described  by  Armendariz  et  a  I .  H8D. 
Briefly,  the  area  is  composed  of  hillocks  2  to  3  meters  in  height  and 
randomly  spaced  approximately  5  to  10  meters  apart.  The  distance  constant 
of  the  wind  instruments  is  1.3  meters  and  data  were  automatica1 ly  col¬ 
lected  at  a  rate  of  one  sample  per  second.  The  thermocouples,  r.0025  cm 
in  diameter  and  made  of  copper  constantan,  have  a  time  constant  less  than 
0.5  seconds  Wind  instruments  were  placed  at  1.5  and  4.0,  or  4.0  and  16.0 
meters  in  neight  on  alternate  poles  in  the  array.  Temperature  differences 
were  recorded  between  0.5  and  4.0  meters  and  between  4.0  and  16.0  meters. 
Data  were  collected  during  the  months  of  January,  February,  and  March  of 
1970.  In  general,  the  collection  of  data  was  made  over  two-hour  periods 
during  both  day  and  night.  There  were  some  periods  of  continuous  data 
collection  lasting  six  or  seven  hours. 


OBSERVED  RElAj  iONSHIPS  BETWEEN  5 !  A !  !C  Si^^LiTY  anju 
SELECTED  DIFFUSION  PARAMETERS 


Relationships  between  meteorological  variables  such  as  static  stability, 
variations  in  wind  direction,  wind  speed,  eic.,  and  atmospheric  diffusion 
parameters  are  complicated  even  in  relatively  simple  situations,  and  cannot 
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T-array  data  presented  here  is  aimed  at  illustrating  some  of  the  complexi¬ 
ties  involved  as  well  as  some  of  the  relationships  between  static  stability 
represented  as  the  difference  in  temperature  between  two  heights  and  certain 
parameters  associated  with  atmospheric  diffusion.  The  objective  is  to 
demonstrate,  for  the  period  of  study  chosen,  the  extent  to  which  a  difference 
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in  temperature  between  two  levels  near  the  ground  and  the  ratio  z/L  are 
related  to  other  parameters  indicative  of  turbulence  and  diffusion. 


Intensity  of  iurbulence 


I.  Longitudinal —  The  longitudinal  intensity  of  turbulence  is 
represented  oy  ou,  which  is  the  standard  deviation  of  the  fluctuation 
along  the  mean  wind  flow.  The  relationships  between  z/L,  aT,  and  °u 
for  selected  heights  are  shown  in  Figure  2.  The  difference  in  temperature 
between  the  indicated  heights  in  each  figure  is  shown  on  the  ordinate, 

°u  on  the  abscissa,  while  z/L  is  the  third  variable  plotted  in  each  figure. 
I  sop  laths  are  drawn  for  values  of  z/L,  and  a  line  was  drawn  by  eye  to 
represent  the  reial ionship  between  AT  and  au. 


In  Figure  2a,  oi  :s  a  function  of  aT  but  the  relationship  is  not  linear. 

It  can  be  seen  f^hat  oM  initially  increases  as  AT  decreases,  but  then  de¬ 
creases  as  AT  decreases  f urther.  The  ratio  z/L  decreases  as  aT  decreases, 
but  the  decrease  is  not  regular  when  z/L  becomes  less  than  approximately 
-0.02.  In  Figure  2a  neither  AT  nor  z/L  is  a  good  indicator  of  ou,  which 
as  we  recall  from  the  theory  presented  in  Section  II  is  directly  related 
to  the  diffusive  power  of  the  atmosphere  for  a  given  mean  wind  speed. 
Figure  2b  differs  from  Figure  2a  only  in  that  z/L  and  ou  are  computed 
at  a  height  of  4  meters  rather  than  1.5  meters.  The  results  are  similar 


and  show  some  indication  of  a  slight  increase  in  ou  with  a  decrease  in 
AT.  As  in  Figure  2a,  there  is  a  poor  relationship  between  z/L  and  ou 
when  z/L  becomes  negative  (less  than  -0.05). 


In  Figure  2c,  where  AT  is  determined  between  4  and  16  meters,  and  z/L 
and  au  at  4  meters,  ou  is  related  to  AT  except  when  aT  is  large  and 
negative.  The  ratio  z/L  decreases  as  aT  decreases  and  ou  increases 
except  at  large  negative  AT.  When  AT  reaches  values  of  -0.5  and  less 
it  no  longer  indicates  the  level  of  turbulence.  In  this  range,  z/L 
increases  as  Oy  increases.  When  AT  exceeds  -0.3,  both  AT  and  z/L  are 
related  to  ou. 

The  fact  that  relationships  between  the  parameters  considered  here  change 
drastically  at  different  levels  in  the  atmosphere  is  illustrated  in  figure 
2d.  This  figure  is  similar  tc  Figure  2c  with  the  only  change  being  that 
z/L  and  au  are  determined  at  16  ineters  rather  than  al  4  meters.  In  figure 
2d,  ou  increases  as  AT  decreases,  but  for  any  At,  ou  may  vary  from  approxi¬ 
mately  1.5  to  3  mps.  The  relationship  between  AT  and  au  is  quite  weak. 
z/L  and  AT  show  a  simi lar  relationship  to  au. 

2.  Latera I  —  In  Figure  3  the  lateral  intensity  of  turbulence,  av, 
is  shown  as  a  function  of  z/L  and  aT  for  the  same  data  as  in  Figure  2. 

The  relationships  between  the  variables  in  Figure  3a  are  the  same  as  those 
discussed  above  for  Figure  2a.  By  comparison,  Figure  3b  differs  con¬ 
siderably  from  its  counterpart,  Figure  2b.  Here,  ov  decreases  wi  it. 
decreasing  aT  when  aT<-0.5.  As  a!  becomes  smaller,  av  varies  by  a  factor 
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Figure  2.  GRAPHICAL  RELATIONSHIPS  BETWEEN 
AT,  (Tu,  AND  z/L. 
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Figure  3.  GRAPHICAL  RELATIONSHIPS  BETWEEN 
AT.  o„,  AND  z/L. 


of  2-5  for  3  given  vs \ us  of  A*  •  Lxcepf  for  nesr— neulrsl  conditions*  z/L 
is  not  a  good  indicator  of  ov  and  when  z/L<-0.l,  the  relationship  between 
z/ •-  and  Cy  apparenriy  nas  vanished* 

In  Figure  3c,  AT  is  not  related  to  oy.  ov  varies  from  approximately 
1.5  to  3.5  for  any  value  of  AT.  The  ratio  z/L  also  is  not  related  to 
oy  any  better  than  to  AT.  These  rather  poor  relationships  extend  to 
Figure  3d,  which  shows  that  ov  is  not  related  to  AT  and  tha+  ov  varies 
between  I  and  4  mps  for  any  given  value  ox  AT.  In  addition,  ther  .  is  no 
apparent  relationship  between  z/L  and  av. 

3.  Vertical  —  Relationships  between  z/L,  AT,  and  ow  are  shown  in 
Figure  4.  Conditions  represented  in  this  figure  are  similar  to  those  in 
Figures  2  and  3.  In  Figure  4a,  cw  is  a  function  of  AT>0.  When  AT  is 
less  than  this  value,  ov  does  not  depend  upon  AT.  The  intensity  of  the 
vertical  component  of  turbulence,  ow,  increases  as  z/L  decreases  to  values 
ne'r  -0.02,  but  beiow  Ihis  value  z/L  and  ow  are  essentially  constant. 

Figure  4b  is  similar  to  Figure  4a  and  shows  that  aw  increases  with  a  de¬ 
crease  in  aT>-0.4,  and  decreases  for  smaller  values  of  AT.  in  this  figure, 
z/L  is  a  poor  indicator  of  ow  for  all  values. 

Figure  4c  is  somewhat  similar  to  Figures  4a  and  b  and  shows  aw  to  be 
related  to  AT,  but  the  relationship  is  not  linear.  ow  increases  as  AT 
decreases  to  values  above  about  -0.3,  but  as  AT  becomes  smaller  ow  remains 
essentially  constant,  with  some  slight  tendency  to  decrease  as  AT  approaches 
-1.0.  The  ratio  z/L  decreases  with  an  increase  in  ow  for  all  AT  except 
for  AT  <-0.3,  in  which  case  larger  negative  values  of  z/L  are  associated 
with  smaller  values  of  aw. 

Figure  4d  differs  from  Figures  4a-c  in  that  ow  increases  as  AT  decreases 
throughout  the  entire  range  of  AT.  In  addition,  z/L  decreases  as  ow  in¬ 
creases  and  AT  decreases.  In  this  figure,  both  AT  and  z/L  are  related  to 
ow  in  a  reasonably  definite  way. 

B.  Variability  of  Wind  Direction 


The  variability  of  wind  direction  is  considered  in  Figure  5.  The  condi¬ 
tions  in  this  figure  are  similar  to  those  in  Figures  2-4  except  that  the 
standard  deviation  ot  component  wind  speeds  has  been  replaced  by  cr^.  As 
we  can  see  from  Eq.  (10),  is  a  measure  of  the  intensity  of  turbulence 
for  a  given  mean  wind  speed.  In  Figure  5a,  a*  is  not  related  to  AT ;  in 
addition,  z/L  is  not  related  to  c?A  in  any  definite  way.  When  z/L  ^-0.02, 
oA  may  vary  by  a  factor  of  about  4  while  z/L  remains  essentially  constant. 
Similar  relationships  are  observed  in  Figure  5b  except  that  when  z/L 
becomes  <-0.l,  larger  values  of  a ^  are  associated  with  smaller  values  of 
z/L. 
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Figure  3r  is  similar  to  Figure  5b-  Again,  o*  is  not  a  function  of  41 
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AT  of  approximately  -0.8C.  For  z/L  <-0.04,  increases  as  z/L  decreases. 
Thus,  AT  is  not  a  good  indicator  of  o^,  while  z/L  is  valid  only  at  values 
<0.04.  Simi lar  relationships  are  observed  in  Figure  5d  except  that  when 
z/L  approaches  -O.i,  a ^  fends  to  increase  as  z/L  decreases. 


Figures  6a  and  6b  shot.  a %  at  1.5  meters  versus  &T  between  0.5  and  4.0 
meters  and  5-5  and  4.0  meters  for  the  period  January-J4arch  1970.  Figure 
6a  differs  considerably  from  Figure  6b  with  the  difference  due  entirely 
to  a  slight  change  in  the  layer  over  which  &T  was  measured.  In  excess 
of  150  hours  of  data  are  plotted  in  each  part  of  these  figures  with  each 
point  representing  a  1 5-ir.i nute  average  value.  This  larger  sample  of  data 
agrees  in  general  with  the  results  from  the  much  smaiier  sample  considered 
in  Figure  5.  The  major  difference  is  that  the  range  of  oa  tor  a  given  AT 
is  much  larger  in  Figures  C?  and  6b  than  in  Figure  5.  Record  et  al.  K) 
found  04  to  vary  with  wine  speed  during  stable  and  unstable  conditions; 
huwSvs r j  this  rc!stion$hip  «2S  oyam » non ^ 


From  the  above  discussion  it  is  evident  that  z/L  and  &T  are  somewhat 
related.  Of  course,  they  should  be  since  AT  appears  in  the  Richardson 
number  (Ri )  and  z/L  =  Ri  for  near-neutral  cond'ticns.  The  experimental 
relationship  is  shown  in  Figure  7  for  the  period  january-iferch  i970.  The 
relationship  is  quite  poor.  The  scatter  when  |z/L|  and  |aT|  are  large 
suggests  that,  at  most,  one  or  possibly  neither  of  these  parameters  can 
be  used  as  an  indicator  of  the  level  of  turbulence,  in  principle,  z/L 
should  be  the  better  indicator;  however,  some  refinements  in  its  compu¬ 
tation  may  be  necessary  when  atmospheric  conditions  vary  greatly  from 
neutra ! . 


C.  Wind  Speed 


The  relationships  between  the  longitudinal  (u)  and  lateral  (v)  component 
wind  speeds  at  1.5  meters,  and  T  oatween  0.5  -  4.0  and  1.5  -  4.0  meters, 
are  shown  in  Figures  8a  and  8b  for  the  period  January-March  1970.  For 
both  wind  components,  the  average  wind  speed  increases  as  AT  decreases 
for  Ai  >  0.  When  aT  <  0,  the  functional  dependence  disappears.  A  line 
has  been  drawn  by  eye  through  the  points  on  each  part  of  the  figure  when 
Al  >  0.  The  range  in  wind  speed  for  a  given  AT  is  large,  particularly 
for  the  lateral  component,  but  there  is  Istt 
present  in  the  data  when  AT  is  positive. 
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CONCLUSIONS 

A_. _ For  the  Layer  1.5  -  4.0  Meters 

The  standard  deviations  of  a! I  three  components  of  the  wind  vector  at 
1.5  and  4.0  meters  (ou,  ov,  ow)  generally  increase  with  a  decrease  in  AT 
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Figure  6(a).  GRAPHICAL  RELATIONSHIPS  BETWEEN 
AT(0,i>-4.0  m)  AND 


Figure  6(b). 


GRAPHICAL  RELATIONSHIPS  BETWEEN 
AT(1. 5-4.0  m)  AND  0/\. 
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Figure  7.  GRAPHICAL  RELATIONSHIP  BETWEEN  AT  AND  z/L. 
(THE  DOUBLE  LINE  INDICATES  A  BREAK  IN  THE  ABSCISSA.) 
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(1-5  -  4-0  meters)  when  AT  >  -0.25C,  but  when  AT  <-0.25C,‘ the  magnitude 
of  the  o's  generally  becomes  smaller  as  AT  becomes  more  negative.  z/L 
also  decreases  as  the  o's  increase  for  values  o-'  z/L  >  -0-02,  but  for 
smaller  values  of  z/L  no  general  relationship  appears  to  hold.  The 
standard  deviation  of  wind  direction  is  net  related  to  AT  or  z/L. 

B.  For  the  Layer  4.0  -  16.0  Meters  ' 

Standard  deviations  of  the  lateral  component  (cv)  and  wind  direction 
(0^)  are  not  related  to  AT  or  z/L  for  any  range  of  values.  The  standard 
deviations  Oi  the  longitudinal  and  vertical  components  (c?u  and  ow)  in¬ 
crease  with  a  decrease  in1  AT  (4.0  -  16.0  meters)  for  AT  >  -0.25G,  but 
when  AT  <  -0.25C,  the  magnitude  of  the  o’s  either  remains  constant  or 
becomes  smaller  as  AT  becomes  more  negative.  At  the  top  of  the  layer 
(16.0  meters),  ou  and  ow  tend  to  increase  with  a  decrease  in  AT  although 
the  relationship  is  poor. 

C.  Wind  Speed  vs  Stability 

i hs  wind  speed  increases  as  The  cegrce  ct  slat:_  sfaL:!it/  decreases, 
but  only  during  stable  conditions.  During  unstable  conditions  (de¬ 
crease  of  temperature  with  height),  the  wind  speed  may  vary  from  1-8  m/sec 
and  is  independent  of  the  degree  of'  instability. 

D.  Relationship  Between  z/L  and  AT 

z/L  increases  as  AT  increases  in  near-neutral  conditions,  but  the 
relationship  diminishes  when  the  degree  of  stability  is  large  cr  Smaii. 

i 

E.  Stability  vs  Diffusion  Pates 


If  we  assume  that  the  rate  of  diffusion  is  determined  by  the  intensity 
of  turbulence  (fluctuations  about  the  mean),  AT  measured  through  a 
shallow  layer  near  the  ground  is  not  a  good  indicator  of  diffusion  rates, 
particularly  when  the  atmosphere  is  statically  unstable  (AT  negative).. 
z/L  does  not  appear  to  be  a  significantly  better  indicator  man  AT  during 
highly  unstable  conditions. 
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